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Abstract 

TAUOLA is a Monte Carlo generator dedicated to generating tau-lepton decays and 
it is used in the analysis of experimental data both at B-factories and LHC. TAUOLA 
is a long-term project that started in the 90’s and has been under development up to 
now. In this note we discuss the status of the predominant hadronic tau-lepton decays 
into two (Br ~ 25.52%) and three pions (Br ~ 18.67%). 

1 Introduction 

The precise experimental data for tau lepton decays collected at B-factories (both BELLE [1J 
and BABAR [2]) provide an opportunity to measure the Standard Model (SM) parameters, 
such as the strong coupling constant, the quark-mixing matrix, the strange quark mass etc, 
and for searching new physics beyond SM. The leptonic decay modes of the tau lepton allow 
to test the universality of the lepton couplings to the gauge bosons. Its hadronic decays (the 
tau lepton, due to its high mass, is the only one that can decay into hadrons) give information 
about the hadronization mechanism and resonance dynamics in the energy region where the 
methods of the perturbative QCD cannot be applied. Also hadronic flavour-violating and 
CP violating decays of tau lepton allow to search for new physics scenarios. 

Hadronic tau lepton decays are also a tool in high-energy physics. In the last two years 
tau physics has received great attention after the discovery of the Higgs boson by the CMS 
and ATLAS Collaborations pm- While this particle resembles in most features the SM 
Higgs boson, the data, showing an enhancement of the two-photon event rates, as well as 
a reduced rate into tau and W pairs, could indicate extended symmetries. Another yet 
unsolved question is related to the origin and precise nature of violation of symmetry under 
charge and parity transformations, the so called CP violation problem. 

Since the 90’s TAUOLA is the main Monte Carlo (MC) generator to simulate tau- 
lepton decays (5]. It has been used by the collaborations ALEPH [6], CLEO [7], at both 
B-factories (BABAR [2] and BELLE CD as well at LHC pm experiments 0 for tau decay 
data analysis. The library TAUOLA can be easily attached to any MC describing the 
production mechanism like KORALB, KORALZ and KKMC [8| [9l HO]. The code provides 
the full topology of the final particles including their spin. 

1 At LHC, at the moment, tau decays are only used for identification and are not used for studying their 
dynamics. However, the dynamics of tau decays are important for both modeling the decays (and therefore 
the reconstruction and identification) and for measuring the decay polarization. 



In view of the forthcoming BELLE-II project mi it is important to revise the TAUOLA 
context in detail and we start onr discussion by reviewing the models describing the tau- 
lepton decays into two (Br ~ 25.52%) and three pions (Br ~ 18.67%). 


2 Two-pion and three-pion hadronic currents 


For r decay channels with two mesons, A {j )\) and h 2 (p 2 ), the hadronic current reads 


J 11 = N[( Pl -p 2 rF v (s ) + (pi +p 2 y i F s (s)], (1) 

where s = (p\ + p 2 ) 2 , F v (s) and F s (s ) are the vector and scalar form factors, respectively. 

In the general case, a hadronic current of a two-meson tau-lepton decay mode depends 
on both vector ( F v ) and scalar ( F s ) form factors. However, in the isospin symmetry limit, 
m r = m n o , the scalar form factor vanishes for the two-pion decay mode and the current is 
described by the vector form factor only. The normalization coefficient (IV) depends on the 
mode and equals 1 for the two-pion mode. 

Currently MC TAUOLA includes four parametrizations for the vector form factor of two 
pions F^(s): 

• Kuhn-Santamaria (KS) parametrization [12] : 
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where M is a resonance mass and T 7rJr (s) is the resonance energy-dependent width 
that takes into account two-pion loops; 


• Gounaris-Sakurai (GS) parametrization used by BELLE [T], ALEPH and CLEO col¬ 
laborations: 
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where f(s) includes the real part of the two-pion loop function; 
parametrization based on the Resonance Chiral Lagrangian (RChL) [13]: 
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where B 22 is the two-meson loop function For the physical meaning of the model 
parameters F v . and Gy i see mm-, 

• combined parametrization (combRChL) that applies dispersion approximation at low 
energy and modified RCliL result at high energy [15] : 
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where So is the high energy limit of the dispersion representation applicability. It is 
supposed to satisfy l.OGeV 2 < so < 1.5GeV 2 [15] and we leave it to be fitted. 


In all the above parametrizations, except for the RCliL one, the pion form factor is given 
by interfering amplitudes from the known isovector meson resonances p(770), p'(1450) and 
p"(1700) with relative strengths 1, /3 and 7 . Although one could expect from the quark 
model that f3 and 7 are real, we allow the parameters to be complex (following the BELLE, 
CLEO and ALEPH analysis) and their phases are left free in the fits. I 11 the case of the 
RCliL parametrization we restrict ourselves to the p(700) and p^HSO) contributions, with 
the relative p' strength (that is a combination of the model parameters F Vi , Gy i and F) 
being a real parameter. 

For the energy-dependent width of p(770)-meson two-pion and two-kaon loop contri¬ 
butions are included for both RChL and combRChl parametrizations through the loop 
functions B 22 an d B 22l whereas in the case of KS and GS the p width is approximated only 
by the two-pion loops through The p'(1450) and p"(1700) widths include only two-pion 
loops for all parametrizations except for the RChL one. In the case of the RChL one both 
two-pion and two-kaon loops are included. 

Results of the fit to the BELLE data [I] are presented in Figs. 1 and 2. The best fit 
is with the GS pion form factor (y 2 = 95.65) and the worst one is with the RChL one 
(y 2 = 156.93), which is not able to reproduce the high energy tail. As it was mentioned 
above, the two main differences of the RChL parametrization compared to the others are 1) 
the p"-meson absence, 2) a real value of the p / (1450)-meson strength. To check the influence 
of the p // (1700) on the RChL result, the p"(1700) resonance has been included in the same 
way as it was done for p'(1450); however, this inclusion has not improved the result (see, 
Fig. 2, the right panel). 


2 Comparing the imaginary part of the loop function B 22 \ 
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In an effective field theory, like RChL, complex values come only from loops. Therefore, 
we conclude that missing loop contributions could be responsible for the disagreement and 
that the complex value of the /3 and 7 parameters might mimic missing multiparticle loop 
contributions. The same conclusion was reached in Ref. |T3] where it was stressed that the 
two-pseudoscalar loops cannot incorporate all the inelasticity needed to describe the data 
and other multiparticle intermediate states can play a role. We are going to check this point 
by adding first a four-pion loop contribution to the p'-resonance propagator. 




Figure 1: The pion form factor fit to Belle data [I]: the GS parametrization (left panel), the 
RChL parametrization (right panel). At the bottom of the figure, the ratio of a theoretical 
prediction to the data is given. 

I 11 spite of the fact that the GS parametrization is able to reproduce the BELLE data, it 
has at least two misleading points. One of them is the complex strengths mentioned above, 
the other is the neglect of the two-kaon loops even for the p-resonance propagator. The last 
point will be the object of the further investigation. 

I 11 the case of the combRChL parametrization the fitting curvature demonstrates not a 
smooth behaviour near s = s 0 . Therefore more sophisticated fitting techniques will have to 
be implemented. 

When the r lepton decays into three hadrons and a neutrino, the predominant decay 
mode involves three pions. In the general case Lorentz invariance determines the decompo¬ 
sition of the hadronic current for a three-hadron final state in terms of five Lorentz invariant 
structures [IBj multiplied by hadronic form factors (iq) 

J 1 * = N {T^[d(p 2 - p 3 yF 1 (q 2 ,s 1 ,s 2 ) + c 2 (p 3 -p 1 )' / F 2 (q 2 ,s 1 ,s 2 ) + c 3 {p 1 - p 2 ) u F 3 (q 2 , s u s 2 )] 

+C 4 <f F 4 (<? 2 , Sl, S 2 ) - 47r 2 F 2 c 5 Y vpaPlPlPs F 5 iY , Si, S 2 )}, (2) 

where as usual T /IL/ = g^ — q^qv/q 2 denotes the transverse projector, and q 11 = (P 1 +P 2 +P 3 Y 
is the total momentum of the hadronic system with four-momenta p 4 , p 2 and p 3 and the 
two-pion invariant mass squared is Sj = (77 +p fc ) 2 . 






































Figure 2: The pion form factor within the comRChL parametrization is fitted to BELLE 
data [Li]: the fit without p"(1700) (left panel) and with it (right panel). At the bottom of 
the figure, the ratio of a theoretical prediction to the data is given. 


The scalar functions Fj(g 2 , si, s 2 ) are the three-pion form factors. The vector form factor 
vanishes for the three-pion modes due to the G-parity conservation F 5 = 0. Among the 
three hadronic form factors F t , i = 1,2,3 which correspond to the axial-vector part of the 
hadronic tensor, only two are independent. The pseudoscalar form factor F 4 is proportional 
to m 2 /g 2 [16], thus it is suppressed with respect to F tl i = 1, 2, 3. 

In TAUOLA the following three-pion form factors are available 

• CPC version [5], that includes only the dominant a\ —» pir mechanism production. 
The form factor is a product of the Breit-Wigner for the a\ and p meson; 

• CLEO parametrization; It is based on the Dalitz plot analysis carried out by the CLEO 
collaboration and includes the following intermediate states: a 4 —» (p; p')tt (where 
(p; p') can be either in S- or in D-wave), cp —> a tt, cp —> / 2 (1270)7t, cp —> / 0 (1370)7t. 
In fact, there are two variants of this parametrization. The former is based on the 
CLEO 7i 0/ n 0 7T~ analysis [7] and applies the same current for the 7r - 7r - 7r + mode. The 
later uses the 7r°7r°7r _ current from [7] and the 7r _ 7r _ 7T + current from the unpublished 
CLEO analysis [T8j. It is interesting to point out that a difference between these 
variants of the CLEO parametrization is related with the scalar and tensor resonance 
contributions. All resonances are modelled by Breit-Wigner functions and the hadronic 
current is a weighted sum of their product E|. The model parameters are resonance 
masses and widths as well as the weights; 

3 This approach was contested in Ref. [T5] where it was demonstrated that the corresponding hadronic 
form factors reproduced the leading-order chiral result and failed to reproduce the next-to-leading-order 


one. 

















• modified RChL parametrization |20] - It was based on the RChl results for the three- 
pion currents mi and an additional scalar resonance contribution. The RCliL current 
is a sum of the chiral contribution corresponding to the direct vertex W~ —> mnr, 
single-resonance contributions, e.g. W~ —* p7r, double-resonance contributions, as 
W~ —> gR —y pn. Only vector and axial-vector are included in the RChL hadronic 
currents. More details about the parametrization can be found in |16j . The scalar reso¬ 
nance contribution was included phenomenologically by requiring the RChL structure 
for the currents and modelling the a-resonance by a Breit-Wigner function. 

The CLEO parametrization for the 7r - 7r - 7r + mode has not yet been fitted to the BABAR 
preliminary data |2j, so we used the old numerical values of the parameters fitted to the 
CLEO data [7j. The fit to the BABAR data will be a task for future work. 



Figure 3: The r~ —y decay invariant mass distribution of the three-pion system 

(left panel) and the two-pion system (central and right panels). The BABAR data [2] are 
represented by the data points, with the results from the modified RChL current as described 
in the text (blue line) and the old tune from CLEO [7) (red-dashed line) overlaid. 

The one-dimensional distributions of the two- and three-pion invariant mass spectra 
calculated on the base of the modified RChL parametrization have been fitted to the BABAR 
preliminary data [2j. It must be pointed out that without the scalar resonance contribution 
the RCliL parametrization provides a slightly better result than the CLEO parametrization 
whereas the scalar resonance inclusion strongly improves the low two-pion mass invariant 
spectrum. Discrepancy between theoretical spectra and experimental data can be explained 
by missing resonances in the model, such as the axial-vector resonance a\ (1600), the scalar 
resonance /o(980) and the tensor resonance /2(1270). Inclusion of these resonances in the 
RCliL framework will be a future task. 

Comparison of the 7r - 7r - 7r + current in the framework of the modified RChL with the 
ChPT result has demonstrated that the scalar resonance contribution has to be corrected 
to reproduce the low energy ChPT limit. The corresponding calculation is in progress. 












3 Conclusion 


In this note we have reviewed the physical frame of the Monte Carlo generator TAUOLA 
for the part related with two- and three-pion decay modes. We have briefly described the 
models used for the hadronic currents as well as the results of the preliminary fits to the 
experimental data. 

In the case of the two-pion modes we have made a fit to the BELLE data for the three 
parametrizations of the pion form factor. We have reproduced the data with the Gounaris- 
Sakurai pion form factor parametrizations while have failed with the RChL parametrization. 
We have concluded that the complex value of the resonance strength, used in the Gounaris- 
Sakurai parametrization, could mimic the missing multiparticle loops. This idea will be 
tested by adding the four-pion loops to the //-resonance propagator that will be the object 
of future study. 

In the case of the three-pion mode TAUOLA has two parametrizations for the hadronic 
current. At present we have fitted the modified Resonance Chiral current, which includes 
the cr-meson contribution in addition to the Resonance Chiral Lagrangian current, to the 
preliminary BABAR data. We have demonstrated that the difference between the model 
and the data is within 5 — 7% at the high-tail of the two- and three-pion spectra. This 
discrepancy is related to the neglect of the resonances, like / 0 (980), a\ (1600) and / 2 (1270), 
in the model. The incorporation of these resonances is in progress. 

The same type of comparison will have to be done for the CLEO n~n~n + current. 

One of the key points of the BELLE-II project mentioned in Introduction is tau physics m 
|2T] . The first physics run is planned in the fall of 2018. Both allowed and forbidden tau 
decay modes in Standard Model will be measured. It should record a 50 times larger data 
sample than the BELLE experiment until 2022. Indeed, the two- and three-pion decay 
modes will be used for measuring resonance parameters on one side and on the other these 
processes constitute an ever present background to the hadronic decay modes, forbidden 
in the Standard Model, for example the lepton flavour violated modes, like r —y Inn [22] . 
All these considerations underline the importance of a precise theoretical description of the 
two- and three-pion decay modes. 
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